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ABSTRACT 

The paper considers the numerical solution of the unsteady, 
three-dimensional, Euler equations to obtain the blade surface pressures of an 
advanced propeller at an angle of attack. The specific configuration 
considered is the SR7L propeller at cruise conditions with a 4.6° inflow angle 
corresponding to the +2° nacelle tilt of the Propeller Test Assessment (PTA) 
flight test condition. The results indicate nearly sinusoidal response of the 
blade loading, with angle of attack. For the first time, detailed variations 
of the chordwise loading as a function of azimuthal angle are presented. It 
is observed that the blade is lightly loaded for part of the revolution and 
shocks appear from hub to about 80 percent radial station for the highly 
loaded portion of the revolution. 

INTRODUCTION 

Advanced propeller tests in wind tunnels and in flight have been 
conducted to understand the aerodynamics and acoustics of advanced designs. 
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Aerodynamic tests of the 9-ft diameter, single rotation SR7L NASA/Haml 1 ton 
standard design were run In the ONERA S1-MA wind tunnel In Modane, France. 
Blade surface steady and unsteady pressures were measured during early 1987 on 
a two blade configuration (The complete eight blade propeller was not 
tested due to drive power limitations.) The propfan test assessment flight 
program (1987-1988) was managed by Lockheed Aeronautical Systems Company under 
a contract to the NASA Lewis Research Center. The 9-ft propfan was flight 
tested on a modified, instrumented Gulfstream GII business jet (Fig. 1). The 
objectives of the PTA program were to evaluate the propfan structural 
Integrity, source noise, cabin noise, flyover noise and enroute noise. Thus, 
an extensive data base has been established for the development and validation 
of propfan aerodynamic and acoustic analyses. 

In the PTA flight test program, the propfan propulsion system was mounted 
on the left wing of the modified GII aircraft. This location was chosen to 
avoid flow field Interference with the main propulsion Gulfstream engines and 
to permit systematic Investigation of the propfan characteristics. 4 A nacelle 
tilt arrangement was employed to vary the Inflow angle to the propfan. The 
variations In Inflow conditions were used to evaluate the effects on cyclic 
stress of the propfan over a wide range of operating conditions. The Inflow 
conditions also play an Important role In the generation of noise. The three 
nacelle tilt angles of the test were -3°, -1° (tilt down) and +2° (tilt up). 
The average Inflow angle Is dependent on the airplane angle of attack, propfan 
upwash angle and nacelle tilt. The complete nacelle tilt variation of 5° (-3° 
to +2 8 ) provided an Inflow angle variation of 5°. 5 

Acoustic and performance data were acquired over a wide range of 
operating conditions. The test points were specifically chosen to determine 
the near field noise characteristics as a function of flight Mach number, 
altitude, propfan tip speed, power and inflow angle. 
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The near field acoustic data for the baseline design-point cruise 
(M = 0.8) case of -1° nacelle tilt have been analyzed In detail. 5 The 
acoustic data measured on a wing boom outboard of the propfan show that the 
first five harmonics have maximum sound pressure levels (SPL's) In the plane 
of the propfan. The SPL's decrease In level with Increase In harmonic order 
and all have similar directional characteristics. The sound pressure levels 
were also predicted using Hanson's frequency domain propeller noise radiation 
theory. 6 The predictions show that the SPL levels are maximum at 0.25 
diameter aft of the propeller plane. The levels are overpredicted by about 3 
dB. 

The effects of nacelle tilt on the fuselage and wing boom SPL's were 
studied both at the cruise and takeoff conditions. Increasing nacelle tilt 
resulted In significant reductions In SPL at the fuselage for all tip speeds. 
At the wing boom, the SPL's were found to increase with Increase in tilt 
angle. On the average an acoustic sensitivity of about 1 dB/degree of nacelle 
tilt was observed. Predictions of the SPL's for the nacelle tilt test cases 
were also done by Hamilton Standard; however, the predicted fuselage and wing 
boom SPL's showed very little sensitivity to nacelle tilt angle. 

Prediction of the effect of inflow angle on the sound pressure levels 
requires accurate computation or measurement of the unsteady blade surface 
pressures. The linear analytical methods for computing the unsteady blade 
loading may not predict the loading accurately enough at higher inflow angles 
due to nonlinear processes Involved. In the present paper, the unsteady blade 
pressures are computed by numerically solving the three-dimensional unsteady 
Euler equations. The only published three-dimensional unsteady Euler solution 
for the propeller at angle of attack is that of Wltfield et al.^ The emphasis 
In that work was on validation of the solution algorithm for different 
configurations. The present paper provides, for the first time, detailed 
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unsteady Euler solutions of the propeller at an angle of attack. The angle of 
attack Induced unsteadiness is Illustrated in the form of chordwlse, spanwise 
and azimuthal pressure distributions. 

SOLUTION OF THE UNSTEADY EULER EQUATIONS 
The unsteady three-dimensional Euler equations governing the inviscid 
flow through a propeller are solved using the solution procedure developed by 
Whitfield. 7 * 8 (The details may be found In the references cited. Here only a 
brief mention of the technique will be made.) The equations in conservative 
differential form are transformed from a Cartesian reference frame to a time 
dependent body fitted curvilinear reference frame. The transformed equations 
are discretized employing a finite volume technique. An approximate Rlemann 
solver Is used for block Interface flux definitions. A Lower-Upper (LU) 
Implicit numerical scheme which possesses apparent unconditional stability Is 
used to solve the discretized equations. The flow field is represented by 
multiblock composite grids to overcome computer core memory limitations. The 
solution at each time step Is updated t»y having only one block In memory while 
the other blocks are stored In solid-state storage devices (SSD). 

THE FLOW CONFIGURATION AND COMPUTATIONAL GRID 
The configuration considered here Is the eight blade 9— ft diameter SR7L 
propeller (Fig. 1), of the PTA program. The Inflow angle chosen for the 
present computations is 4.6® (corresponding to nacelle tilt of +2° in the PTA 
program). The direction of rotation of the propeller and the reference for 
the coordinates are shown on the side of the figure. 

Each blade passage is discretized employing a 71 by 33 by 11 grid, in the 
axial, radial and circumferential directions, respectively. Each passage grid 
is divided into two blocks for computational convenience as mentioned before. 
Thus, 16 blocks of grids are employed to describe the entire propeller flow 
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field with 206,184 nodal points. Each blade surface has 35 by 21 (chordwlse 
by spanwlse) grid points with higher resolution near the leading and trailing 
edges. The solution presented here Is for the cruise condition, Mach number = 
0.801 and advance ratio 3.122, with an Inflow angle of 4.6°. 

RESULTS AND DISCUSSION 

The unsteady three-dimensional Euler equations have been solved (from an 
Impulse start) for three complete revolutions of the propeller to get a 
reasonably accurate solution. The results of the third cycle are fairly 
stabilized as can be seen from Fig. 2. The figure shows the variation of the 
total power coefficient with azimuth angle. The total power coe f ficient 
becomes nearly constant during the entire third revolution of the propeller. 
The results of this cycle are analyzed further and presented. The predicted 
total power Is about 8 percent higher than the measured value. No chordwlse 
pressure measurements are available for a direct comparison of the 
predictions. However, the predictions of unsteady surface pressures employing 
the present solution technique were compared with data for SR3 eight bladed 
model propeller in Ref. 7 (Fig. 5). The comparisons were done at blade chord 
location x/c = 0.1 and 0.5 at blade span locations 0.75 and 0.88 radius for 
an inflow angle of 4° at Mach numbers 0.6 and 0.8. At the axial station 
x/c = 0.5 the agreement of the predictions with measurements were very good on 
both suction and pressure surfaces at Mach numbers 0.6 and 0.8. At the axial 
location x/c =0.1, predicted surface pressures showed reasonable agreement 
with data for M = 0.6 on both, suction and pressure surfaces. At the higher 
Mach number, 0.8, the pressure side was in reasonable agreement with the 
experiment but the suction side was not. The discrepancy was attributed to 
the coarse grid near the leading edge for the rapidly expanding flow on the 
suction surface at this Mach number. In Ref. 7, 21 (axial) by 12 (radial) 
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grid points on the blade were employed. With 35 (axial) by 21 (radial) grid 
points on the blade surface and higher resolution near the leading and 
trailing edges, employed In the present study a higher accuracy of the 
predictions Is expected. The grid size and distribution were arrived at 
during the comparison study of steady Euler solutions with measurements for 
the SR7L two bladed configurations over a range of Mach numbers. 9 

Figure 3(a) shows the variation of the single blade power coefficient 
with azimuth angle for two blades, one starting from <|» = 0 (or <|> s = 0) 
position and other starting from <j> = ir (or 4> s = ir). The expected sinusoidal 
variation of the loading due to angle of attack, (a = 4.6°) Is clearly 
exhibited. However, the variation of power during a cycle Is significantly 
higher than one might expect. For the case of 4.6° Inflow angle, the power 
coefficient per blade varies +81 percent during a cycle. A careful look at 
the plot reveals additional details: the loading Is not exactly sinusoidal; 

the loading Is not minimum at <j> - 0 and the maximum In not at <J> * ir as 
originally believed. The rising portion of the wave lasts for about 9/8 ir 
while the falling portion of the wave last about 7/8 ir. For the blade 
starting at $ - 0, the minimum loading occurs at ir/8 and the maximum loading 
occurs at about (ir/8 +9/8 ir = 5/4 ir) or 225°. These points have to be borne 
In mind during the analysis of the results. 

A Fourier transform on the third cycle loading gives the blade power 
coefficient Cp as Cp = 0.238 - 0.1987 cos wt - 0.0511 sin wt. The loading 
spectrum shown In Fig. 3(b) indicates that the first harmonic dominates the 
loading. The oscillating part of the Cp lags behind the motion of the blade 
by 14.4°. 

The reduced frequency, which is a measure of the unsteadiness of the flow 
is defined as, k * wc/2U<x>, where c Is the blade chord length, w is the 


6 


blade rotational frequency and lb Is the free stream velocity. The reduced 
frequency for the present calculations is k = 0.088, based on the chord 
length at 75 percent radius. 

For comparison purposes a steady state solution was obtained for the same 
advance ratio (3.122) and Mach number (0.801) but with zero angle of attack. 
The steady blade surface pressure distributions are shown in Fig. 4. The 
figure shows chordwlse pressure distributions at eight spanwlse locations. 

The unsteady blade surface pressure distributions presented in the succeeding 
figures will be at one or more of these spanwlse locations. Thus, a direct 
comparison of the steady and unsteady pressure distributions can be made. 

The detailed azimuthal variations of the chordwlse loading and blade 
surface pressure distributions are Illustrated In Figs. 5 to 10 for three 
radial stations namely, r/R = 0.36, 0.66 and 0.96. Figure 5 presents the 
chordwlse loading variations as a function of azimuth angle, at the radial 
station r/R = 0.36. At this radial station, the blade Is negatively loaded 
for the front half (0 to 50 percent) of the cord length for about 3/8 of a 
revolution (<j> * 0°, 45°, 90°) and the loading Is maximum at <|> = 225°. Thus, 
the chordwlse loading undergoes a cyclic variation of substantial amplitude. 
The shape of the loading curve also changes enormously during the cycle. The 
reasons for this change become apparent from the actual blade surface 
pressures shown In Fig. 6 for the same radial azimuthal locations. The 
trailing edge shock appears on the suction surface during half the cycle 
(<|> = 135° to 270°) while the shock appears on (or extends to) the pressure 
surface for <j> = 135° to 225°. The dramatic changes in the loading pattern 
are clearly illustrated In this figure. 

Figure 7 shows the azimuthal variation of chordwise loading distribution 
at r/R = 0.66. The part of the revolution during which significant negative 
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loading exists. Is reduced for this radial station. A continuous change of 
the space of the loading curve occurs over the revolution. The second half 
(50 to 100 percent chord) of the blade Is highly loaded from <j> - 90° to 
225°. The blade surface pressures shown In Fig. 8 Indicate the appearance of 
the trailing edge shock for 4> = 90° to 270° and a pressure jump on the 
pressure side occurs only for <J> - 90° to 225°, Thus, for the two radial 
stations (0.36 and 0.66), the shocks appear only during parts of the cycle. 

It Is also Interesting to note that the shock on the pressure surface moves 
upstream as we go from <j> = 90° to 225°. 

The periodic motion of shock waves on oscillating airfoils In transonic 
flow have been extensively studied and documented.^® Three types of shock 
motion have been Identified; (1) Sinusoidal shock-wave motion (type A): The 

shock-wave moves almost sinusoidally and exists during the complete cycle. 

(2) Interrupted shock-wave motion (type B): This motion Is like type A but 

the shock disappears during part of a cycle and (3) Upstream propagated 
shock-wave (type c): At slightly super critical Mach Numbers, periodically a 

shock-wave Is formed on the upper surface of the airfoil and moves upstream 
while Increasing In strength. No such detailed experimental study exists for 
the three-dimensional transonic flows of the type considered in this paper. 

But, the results presented In Figs. 6 and 8 suggest that shock motions similar 
to type B and type C occur in the unsteady three-dimensional transonic flows. 

The azimuthal variation of chordwlse loading at the radial station 
r/R =. 0.96 is shown in Fig. 9. At this radial station, the changes in the 
shape of the loading curve are marginal, compared to those at r/R = 0.36 and 
0.66. However, the location of the maximum loading moves from a point close 
to the trailing edge to a point near the leading edge as we go from 4> = 45° to 
<j> = 315°. The blade surface pressure variations shown in Fig. 10, indicate 


8 


that the shocks have diffused before they reached this radial station. It Is 
Interesting to note that the tip region Is significantly loaded through the 
complete revolution of the blade. 

Figures 11 to 14 show the radial variation of the chordwlse pressure 
distribution for four azimuthal angles, 4> = 0°, 90°, 180° and 270° 
respectively. The radial variations of the chordwlse loading distribution at 
«j> - 0° are shown In Fig. 11. At this azimuthal location the blade Is very 
lightly loaded throughout the span. A small negative loading exists almost up 
to r/R = 0.80, which Is due to the low local flow Incidence angle produced by 
the angle of attack of the propeller. Figure 12 shows the radial variations 
of the chordwise pressure distributions at <|> = 90°. At this azimuthal 
location, shocks exist both on the suction and pressure surfaces from the hub 
region to almost 90 percent span. Up to 80 percent span the loading Is 
primarily on the rearward half (50 to 100 percent chord) of the blade. 

The spanwlse variation of the loading at $ = 180° Is shown In Fig. 13. 

At this location the blade loading experiences a near maximum. At this 
azimuthal location also, shocks appear on both suction and pressure surfaces 
from the hub region. However, on the suction surface the trailing edge shock 
Is diffused beyond about 62 percent span. The shock on the pressure surface 
persists to 80 percent span, as for the <j> = 90° position. Figure 14 shows 
the spanwise variations of the loading for <|> = 270°. At this location no 
shock exists on the pressure side of the blade. On the suction surface, a 
shock appears at about 70 percent chord at the radial station near the hub, 
r/R = 0.36. With increasing spanwlse distance, the shock moves toward the 
trailing edge. The four figures (Figs. 11 to 14) demonstrate that the blade 
loading changes significantly in the chordwlse, radial and azimuthal 
directions during a revolution, due to the angular Inflow to the propeller. 
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The variations of the elemental power coefficient (dCp/dx) for the four 
azimuthal locations (<t> = 0°, 90° , 180° and 270°) are shown In Fig. 15. The 
general shape of the curve is the same as that obtained for a steady flow. 

The magnitude of dCp/dx at any radial distance (x ■ r/R) depends on the 
azimuthal location. Also shown In the figure is the curve for the steady flow 
(a = 0) solution. It is seen that the magnitude of the cyclic variations of 
loading depends on the spanwise location. 

The azimuthal variations of blade pressure coefficient at specified 
points along the chord at the radial station r/R = 0.66 are shown In Fig. 16. 
It can be seen that the response at the Individual points X/C = 0.036, 0.125, 
0.53, 0.63 and 0.84 Is not exactly sinusoidal. The extent of the deviation 
from the sinusoid depends on the chordwlse location and the surface (suction 
or pressure surface) on which the point Is located. The departure from the 
sinusoidal behavior may be due to the high Inflow angle, and the presence of 
shock waves. It should also be noted that the Instantaneous pressure 
coefficients are plotted In this figure rather than the fluctuations about the 
mean value of the pressure. 

CONCLUDING REMARKS 

A numerical solution of the unsteady three-dimensional Euler equations 
has been obtained for the flow through a propeller at an angle of attack to 
the mean flow. The results show the highly nonlinear nature of the loading 
variations during a revolution. For parts of the revolution the blade is very 
lightly loaded, while for other parts It is highly loaded with shocks 
appearing on both surfaces of the blade. The sound pressure levels, the 
directivity and the azimuthal variation will be studied next employing the 
unsteady blade pressure distribution obtained here. 
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FIG. 1 SR7L PROPELLER INSTALLED ON THE TESTBED AIRCRAFT. 
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FIG. 8 AZIMUTHAL VARIATION OF BLADE SURFACE PRESSURES 
AT f/R = 0 . 66 . 
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FIG. 13 SPANWISE VARIATION OF BLADE SURFACE PRESSURES 
AT <P = 180°. 
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